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ABSTRACT 

We use analytical and N-body methods to examine the survival of wide stellar binaries 
against repeated encounters with dark substructures orbiting in the dark matter haloes 
of dwarf spheroidal galaxies (dSphs). Our models adopt cosmologically- motivated con- 
ditions wherein dSphs are dark-matter dominated systems that form hierarchically and 
orbit about a host galaxy. Our analytical estimates show that wide binaries are dis- 
rupted at a rate that is proportional to the local density of dark substructures averaged 
over the life-time of the binary population. The fact that external tides can efficiently 
strip dark substructures from the outskirts of dSphs implies that the present number 
and distribution of binaries is strongly coupled with the mass evolution of individ- 
ual galaxies. Yet we show that for the range of dynamical masses and Galactocentric 
distances spanned by Milky Way dSphs, a truncation in the separation function at 
O'max^O.lpc is cxpected in all these galaxies. An exception may be the Sagittarius 
dSph, which has lost most of is dark matter envelope to tides and is close to full dis- 
ruption. Our simulations indicate that at separations larger than Cmax the perturbed 
binary distribution scales as dN /da ex a~^'^ independently of the mass and density of 
substructures. These results may be used to determine whether the binary separation 
function found in dwarf galaxies is compatible with the scale-free hierarchical picture 
that envisions the existence of dark substructures in all galactic haloes. We show that 
the ACS camera on board of the Hubble telescope may be able to test this prediction 
in dSphs at heliocentric distances ^ 100 kpc, even if the binary fraction amounts only 
10% of the stellar population. 
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1 INTRODUCTION 

Perhaps one of the most striking predictions from the 
present cosmological paradigm (Cold Dark Matter, CDM) 
is the existence of dark matter haloes devoid of baryonic 
material, the so-called "dark galaxies" . Such systems arise 
because primordial gas is unable to cool efficiently by col- 
lisional exicitation of H2 molecules (Haiman et al. 2000), 
whose formation is limited to haloes with virial temper- 
atures above Iff^K (i.e. virial masses ^lO^M©). Haloes 
with masses below ~ 1O*M0 are thus expected to be dark. 
The fact that the CDM halo mass function diverges as 
dN/dm ex m^^'^ (Diemand et al. 2008; Springel et al. 2008 
and references therein) implies that dark galaxies dominate 
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in number over luminous ones. Indeed, we expect of the order 
of 10^^ subhaloes with masses within 1O~^-1O*M0 lingering 
in a Milky Way-size halo (Diemand et al. 2005). 

How can we possibly test the existence of dark galax- 
ies? If dark matter is in the form of a neutral relic particle 
once in thermal equilibrium in the early universe, these ob- 
jects may be lit up by pair annihilation that would result in 
gamma-ray emmissions. Unfortunately, testing this scenario 
on the mass scale of dark galaxies is well beyond the present 
technological capabilities for the vast majority of particle 
candidates (e.g. Diemand et al. 2007). The second alterna- 
tive is to detect dark haloes via gravitational perturbations 
of visible objects. 

In this regard very wide binaries, i.e. those with separa- 
tions ag^ 100 AU, may enable simple experiments to test for 
the presence of dark substructures, for even extremely weak 
tidal perturbations can disrupt them. In the stellar halo of 
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the Milky Way (hereafter MW), perturbers can be identified 
as inhomogeneities in the Galactic potential. Some examples 
are stellar clusters, molecular clouds and compact objects 
(i.e. MACHOS), which altogether cover a mass range that 
overlaps with that of dark galaxies (Carr & Sakellariadou 
1999, Chaname & Gould 2004). Hence isolating the effects 
of different pertubers on binaries is not straightforward. 

Much more interesting targets are binaries in dwarf 
galaxies. The reasons are several: (i) dwarf spheroidal galax- 
ies (hereafter dSphs) have the largest dark matter content 
of any known galaxy type in the Universe (e.g. Mateo 1998, 
Gilmore et al. 2007) . In dSphs the fraction of dark to bary- 
onic mass is so large (10^-10^) that the presence of bary- 
onic inhomogeneities in the dwarf potential can be safely 
neglected, (ii) As a result of their hierarchical formation, 
dwarf galaxies orbiting the MW are expected to contain 
bound substructures. These objects are typically referred 
to as (sub-)subhaloes (hereafter sub^haloes) in cosmological 
simulations. Given that the estimated virial mass of the MW 
dSphs is of the order of Mvir ~ l(f Mq (e.g. Penarrubia et 
al. 2008a, Walker et al. 2009) and that in CDM the halo 
mass function is scale-free, all dark matter substructures or- 
biting in dwarfs are expected to fall below the star formation 
threshold and thus be dark, (iii) Dwarf galaxies are embed- 
ded in haloes with estimated dark matter densities that are 
several orders of magnitude higher than in the solar neigh- 
bourhood. In a scale-free cosmology this implies a higher 
number density of perturbers and thus a higher chance to 
detect their effects on a pre-existing binary population, (iv) 
The stellar populations of dSphs are generally dominated by 
old and intermediate-age components (e.g. Mateo 1998 and 
references therein), which enhances the cumulative effect of 
gravitational perturbations. 

Unfortunaltey, although binaries with intermediate 
( ^ 100 AU) separations have been studied in stellar clus- 
ters with some detail in the past (e.g. White & Ghez 2001; 
Patience et al. 2002), little is known about the number and 
distribution of binary stars in dwarf spheroidal galaxies (e.g. 
Wyse et al. 2002 estimate that 5% of the stars in the Ursa 
Minor dwarf correspond to unresolved binaries from the 
number of blue stragglers detected in this galaxy) . 

This paper is aimed at studying the disruption of wide 
( ^ 100 AU) binaries as a result of interactions with dark 
substructures orbiting in dwarf galaxies, as well as explor- 
ing whether the effects can be detected with our present 
observational capabilities. The arrangement of the paper is 
as follows: §2 outlines the analytical analysis, §3 tests its 
consistency through N-body simulations, §4 applies the re- 
sult to the Milky Way dwarf galaxy population, §5 estimates 
the feasibility of detecting wide binaries in dwarfs and §6 
provides a summary of the results. 



CATASTROPHIC VS. 
ENCOUNTERS 



DIFFUSIVE 



Binary stellar systems travelling through a homogeneous sea 
of perturbers with mass Mp, number density Up and density 
Pp = MpUp suffer the cumulative effect of repeated encoun- 
ters over their life times. Assuming that the relative velocity 
of encounters T/roi is isotropic and follows a Maxwellian dis- 
tribution, the change of energy per unit mass averaged over 



a suite of binary orbits with different orbital phases and 
eccentricities is (Binney & Tremaine 2008; eq. 8.48) 



<A£>=^^^t/(bA.); 



(1) 



where h is the impact parameter and r^ the half-mass radius 
of the perturber. For extended bodies that follow a Hern- 
quist (1990) mass profile U{x) « x^'^ /{I + x)^'^, the point- 
mass approximation, [/ ~ 1, is recovered at a; = h/r^ <L 1. 
For simplicity we shall hereafter assume that the point-mass 
approximation can be applied to all perturbers, whilst leav- 
ing a more stringent analysis for Appendix \K\ 

During the life time of a binary, < A_B > can grow 
(i) as a result of repeated encounters that are not fatal for 
the system or (ii) through single encounters that completely 
shatter the system. In the former case the energy change 
gained after single interactions is small, < /S.E > /\E\ <^ 1, 
where E = —GMi,/{2a) is the binding energy of the binary. 
This is the so-called diffusive regime. Perturbations occur in 
a catastrophic regime if < AE > /\E\ ^ 1. For point-mass 
perturbers the condition for catastrophic encounters is 

1/4 



6< 1.5 



GM;a'' 



(2) 



so that when a tidal perturbation is sufficiently strong to 
be fatal for a binary even at 6,nin/a ^ 1, we consider that 
the effects of single, catastrophic interactions dominate over 
diffusive ones. 

To identify the dominant regime of perturbations we 
define the quantity Afp^min as the perturber mass that yields 
femin/a = 1 in eq. (|2|. Subsequently we scale this equation 
to parameters typically found in dwarf galaxies. Since stars 
in these systems show no rotation and a typical velocity 
dispersion a^, ~ 10 km/s (Mateo 1998; Walker et al 2007) 
the relative velocity of encounters can be assumed to be of 
the order of Vrd ~ V^cr*. Using Mi, = IMq and a = 0.1 pc 
as fiducial binary parameters we find 



Vr, 



( Mi 



1/2 



^/2■ lOkms-i Vl^© O.lpc/ 

so that the boundary between diffusive and catastrophic en- 
counters is placed at Mp ~ IO'^Mq for the range of separa- 
tions relevant in this study. 

In this work we identify "perturbers" with dark mat- 
ter sub^ haloes orbiting in the dark matter haloes of dwarf 
spheroidal galaxies. Recent estimates show that the virial 
masses of the brightest MW dwarfs span a relatively small 
range of values around Mvir ~ lO^M© (e.g. Penarrubia et 
al. 2008a). Adopting the subhalo mass function found in 
CDM simulations, sub^haloes are expected to have masses 
Mp <0.1Afvir « 10® M© > A'/p,niin (Navarro, Frenk & White 
1997; hereafter NFW). Since the expected perturber masses 
are many orders of magnitude larger than Mp.min the effect 
of diffusive interactions can be safely neglected in our study 

Indeed, the fact that interactions between dark matter 
substructures and wide binaries can be treated as catas- 
trophic perturbations is a remarkable result, because in this 
regime the disruption time scales as (e.g. Binney & Tremaine 
2008; eq. 8.51) 



Gpp V 



fGMb\ 



1/2 






(4) 
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Figure 1. Binary distribution as a function of semi-major axis. 
A^i, = 3 • 10^ binaries are generated following an arbitrarily cho- 
sen flat distribution represented as filled circles. Binary and per- 
turber orbits are drawn from an NFW halo N-body model with 
Mvir = IO^Mq and r^ = f.lf kpc evolved during 14 Gyr in iso- 
lation (see text). All binary particles move on orbits with apoc- 
entres smaller than Rf,. The perturber mass spectrum scales as 
dN(Mp)/dMp (X Mp ' , and the total mass in perturbers is 



/sub 

N, 



E, 



Mp^i/Mvir- For these particular models we use 



p — 10^, although these results are independent of our choice 
of Np and Ni,. Arrows show the location of Omax as derived from 
eq. ^. Dashed lines show that for a 2> flmax the separation dis- 
tribution scales dNj,/da oc a~^'^. 



where fccat ^ 0.07 (Bahcall, Hut & Tremaine 1985). Note 
that the disruption of wide binaries is independent of both 
the relative velocity of the encounter and the pertuber mass. 
Thus, observed binary distributions in dwarf galaxies can 
yield significant constraints on the amount and distribution 
of dark matter substructures, pp, in these systems. 
In what follows, it is useful to define the quantity 

{GMt)'^'; (5) 



_ f fccat ^ ' 

"'^''''-{GppAtJ 

so that binaries with initial separations a ^ Omax and a mass 
Mb are expected to be disrupted over a time-scale At. 



3 MONTE-CARLO SIMULATIONS 

To test the consistency of the analytical predictions we simu- 
late the effects of encounters between dark matter subhaloes 
and binaries using a standard Monte-Carlo technique. We 
adopt an "unperturbed" initial separation distribution of bi- 
naries that follows Opik's (1924) law dN^/da oc a^^ within 
lO^-lO'^ AU or, equivalently, that is uniform over the inter- 
val 1.0 < logiQ(a/AU) < 7. For simplicity, we assume U — 1 
and a fiducial binary mass Mt = IMq in eq. ^ (see Ap- 
pendix [X] for simulations where perturbers are treated as 
extended objects). 



The orbital distributions of "stars" and "perturbers" 
are drawn from the N-body models of Peiiarrubia et al. 
(2008b) evolved in isolation for a Hubble time, so that 
At = 14 Gyr in eq. ([5|. These models follow an NFW pro- 
file with a virial mass and concentration Mvir = 1O^M0 and 
Cvir ~ 23.1, respectively. Using the concordance cosmological 
parameters (e.g. Spergel et al. 2007) this choice corresponds 
to scale and virial radii r, — 1.11 kpc and r^ii — 25.6 kpc, 
respectively. 

From the N-body models, Np perturber particles are 
randomly chosen. We also select A'^b binary particles with 
apocentric distances to the dwarf centre smaller than Rb, 
which accounts for the fact that stars are spatially segre- 
gated within the dark matter haloes of dSphs (Penarrubia 
et al. 2008a). Given that the typical half-light radii of bright 
(Mv < -8) dSphs is ~ 300 pc (Mateo 1998), we choose 
a fiducial Rb = 600 pc, i.e. Rb/vs = 0.54. Subsequently, 
we follow for a Hubble time the relative position between 



the star "i" and the perturbers, Vij 



Vb,t 



'p,ji 



where 



l,...,7Vp and i = l,...,A'^i,. Whenever the quantity 
|ri^j| = b finds a minimum we store the parameters of the 
encounter {h,Vrc\)i,j- 

Perturber masses are calculated adopting the subhalo 
mass function typically found in CDM N-body simulations, 
which scales as dN{Mp)/dMp oc M~", where a Ki 1.9 
(Springel et al. 2008). The fact that a < 2.0 is crucial be- 
cause, although the number of perturbers diverges at low 
mass, the total mass in perturbers is dominated by a few 
massive substructures with Mpi^O.lMvir- Here we derive 
the normalization of the perturber mass spectrum by defin- 
ing the quantity fsuh = Z]i=i,...,iVp Mp,i/Mvir. 

Fig. [l] shows that dark matter substructures efficiently 
destroy binaries with long semi-major axes, and that this 
process clearly strengthens as the mass in dark substruc- 
tures, indicated by /sub, increases. The result is a final sep- 
aration distribution that is truncated with respect to the 
initial profile. Remarkably, the perturbed part of the distri- 
bution can be accurately approximated by dNb/da oc a~^'^ 
for ajjamax, independently of perturber mass and densit-al 

For each simulation we mark the value of Omax derived 
from eq. ([5| with arrows. We define the (local) density of 
perturbers as 



(pp) = 



/sub-^'^vir 

Atv/3RI 



(6) 



Interestingly, the analytical value of ttmax accurately sig- 
nals the binary separation at which deviations from the 
unperturbed distribution start to become obvious (see also 
Fig.lAll. 



In the following Section we apply these results to Milky 
Way dSphs in order to examine whether their binary popula- 
tions can probe the existence of dark substructures orbiting 
within their dark matter haloes. 



^ Binary populations experiencing perturbations in a diffusive 
regime also show power-law distributions. However, the slope is 
not unique and depends on several encounter parameters (Yoo, 
Chaname & Gould 2004) 
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Figure 2. Evolution of the sublialo mass fraction /sub = 
y^ Mgub/Mvir measured at different distances from the centre of a 
Milky Way-size halo, i.e. Mvir = W^^Mq (realization taken from 
Pefiarrubia et al. 2010). The infall of massive substructures can 
be easily identified as jumps in the value of /aub- Although this 
quantity suffers strong variations, the time-average values (thick 
lines) remain fairly constant in time. 



Figure 3. Expected truncation in the separation function Omax 
of a binary population with A/;, = IMq that was formed at t^ = 2 
Gyr after the Big Bang (see text). Upper and iotuer panels assume 
that dark structures can be treated as point-mass and extended 
(see Appendix [A]| perturbers, respectively. The horizonthal dot- 
ted lines marks a separation of 1 pc, which roughly corresponds 
to the largest separation expected in a binary sample. 



4 BINARY DISRUPTION IN MILKY WAY 
DSPHS 

Estimating the time-averaged density of substructures in 
dwarf galaxies orbiting about a host is not exempt from 
complexities. In field haloes the process of accretion, tidal 
stripping and merger of small substructures proceeds un- 
interruptedly in a way that keeps the time-averaged local 
subhalo density practically constant in time This is illus- 
trated in Fig. [2J where we show the evolution of the subhalo 
mass fraction /sub(< r) = "^ Msub{< r)/M^ii calculated at 
different radii from a Milky Way-size halo realizatiorrl Infall 
of massive substructures can be easily identified as "jumps" 
in the value of /sub(< ^vir) that propagate towards smaller 
radii as massive substructures sink into the inner-most re- 
gions of the host through dynamical friction (thin lines). 
Remarkably, although the local mass in substructures suf- 
fers strong variations, the time-averaged value of /sub (thick 
lines) stays relatively constant at all radii. 

However, this process is interrupted when field galaxies 
are accreted into a larger host and become satellites: On the 
one hand, the accretion of external sub'^haloes stops after 
satellite galajcies cross the virial radius of the parent galaxy. 
Hence, substructures that are tidally disrupted stop being 
replenished. On the other hand, a fraction of the existing 
substructures is lost after each pericentric passage through 
tidal stripping. The combination of both effects tend to lower 
the mass fraction of sub^haloes in satellites in a monotonic 



^ Since in ACDM the hierarchical formation of haloes is scale- 
free, we expect a similar evolution also on the mass scale of dSphs. 



fashion and induce a correlation between the present value 
of /sub and the location of satellites within the main halo. 

Dealing with these aspects of the problem in a consis- 
tent way requires the aid of cosmological N-body simula- 
tions with a spatial and time resolution that, unfortunately, 
go beyond the present numerical capabilities. Here we make 
three crude assumptions that allow us to estimate the loca- 
tion of the truncation in the binary separation function of 
MW Sphs in a simple way. First we use eq. ([6| to estimate 
the time-averaged local density of perturbers. Second, since 
the amount of stripping experienced by the MW dSphs is 
difficult to gauge from observational data (Pefiarrubia et al. 
2009), we use the statistical correlation between /sub and 
Galactocentric distance found in the Aquarius run (Springel 
et al. 2008) to guide our estimates^ which can be expressed 
as fsuh{R) = fo{R/RvirY for R < -Rvir, where R is the 
Galactocentric distance of the dwarf, /o — 0.2 and u ~ 2.0, 
with a log normal scatter of AlogigC" = 0.5 at fixed radii. 
Third, we adopt Walker et al. (2009) results, who show that 
the estimated dynamical masses of all MW dSphs are com- 
patible with these systems being embedded in dark matter 
haloes with Mvir ~ 10^ M©. 

Under these assumptions, the expected value of amax 
for binaries with Mb = IMq is shown in Fig. [S] assum- 
ing that dark substructures can be treated as point-mass 
{upper panel) and extended {lower panel, see Appendix [A|) 



^ Note that the Aquarius models do not include baryons. Since 
discs strongly enhance the tidal stripping of satellite galaxies 
(D'Onghia et al. 2010, Pefiarrubia et al. 2010), this correlation 
may be stronger in spiral galaxies. 
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separation binary populations in Milky Way dSphs are fea- 
sible with existing instrumentation. 

Let the total number of stars in an observed part of the 
sky be 



log^f^iHeliocentric Distance flkpc) 

Figure 4. Upper panel: Effective signal to noise ratio of bina- 
ries with respect to random pairs S/N = N^,ii^/\/N^,i^ -(- Nmn 
in dSphs with different surface brightnesses (/i) as a function of 
heliocentric distance. We adopt the angular resolution and (a sin- 
gle) field of view of the ACS camera on board of the HST. Lower 
panel: Interval of separations (ai,a2) that maximizes the S/N 
ratio. Note that ai is simply the telescope resolution times helio- 
centric distance. 



perturbers acting on the binary population. Error bars ac- 
count for the scatter found in the cosnrological relation be- 
tween /sub and Galactocentric distance. For ease of refer- 
ence, we add a horizonthal dotted line to mark a separation 
of 1 pc, which roughly corresponds to the maximum sepa- 
ration found in (spectroscopically confirmed) binaries in the 
stellar halo of the MW (Quinn et al. 2009). 

This Figure clearly shows that all dwarf galaxies in the 
MW are expected to have truncated binary separation func- 
tions within the context of cosmologically-motivated CDM 
haloes. Focusing on the extended perturber models, trun- 
cation may be visible at separations as small as ~ 10~^ pc 
for most of the ultra-faint dSphs (Belokurov et al. 2007, 
2009), which have stellar sizes Rb = 2 x rhait ^ lOOpc. 
"Classical" dwarfs on the other hand are more extended, 
500;$ (-Rb/pc) ^800, and tend to reside at further Galac- 
tocentric distances (^80 kpc). For these systems our es- 
timates suggest that the binary population should have a 
separation function truncated beyond ;$0.1 pc. A notable 
exception however is the Sagittarius (Sgr) dwarf, which is 
currently close to full tidal disruption (Niederste-Ostholt et 
al. 2010). Putting its Galactocentric distance in the Aquar- 
ius relation yields a very low content of substructures, /sub — 
2 • 10~*, which likely implies an unperturbed (flmax ~ 1 pc) 
binary function. 



5 OBSERVATIONAL PROSPECTS 

The analysis of the preceding Section suggests that the sep- 
aration of binaries in dSphs may scale as dN/da oc a~^'^ for 
a^O.lpc. Here we examine whether observations of large 



N^ ^flfiNt, 



(7) 



where Q is the survey area, fi is the surface brightness of the 
dwarf galaxy and Nt = J^ f{M) dM is the total number 
of stars per solar luminosity in a dwarf galaxy with distance 
modulus (m — Af ) and stellar luminosity function /(M). 

The number of stars which are actually detected (i.e. 
above the detection limit) is 



Nobs = niiNtg{Am), 



(8) 



where (;[Am] — J_ '" f{M) dM/Nt, with Am = mum — {rn — 
M), denotes the fraction of stars above the limiting magni- 
tude miim. 

Thus the number of physical binaries observable in a 
survey (A^bin) can be written as 

A^bin = giAm)Nobs / -P(a) da = Q fi Ntg''^[Am] / P{a) da; 

(9) 
where P{a) da is the fraction of binaries with separations 
in the interval a and a + da; a^ is the minimal separation 
we are able to resolve in our survey and 02 is the maximal 
separation. 

The observability of wide binaries depends on the con- 
tamination by random stellar pairs. Assunring that the den- 
sity of the stars in the dwarf galaxy is significantly higher 
than the density of foreground/background object, the num- 
ber of random pairs with the range of separations ai < a < 
02 can be easily estimated from the surface density of ob- 
served stars Nobs/^ and the angular area enclosed within 
the annuli ai/D and a2/D 



A'^ran ~ TvCl^ Nt Q [Am 



{al 



al) 



D2 



(10) 



where D is the distance to the galaxy. 

The detectability of binaries in a given separation range 
is determined by an effective signal to noise ratio S/N — 
-'Vbin/V-^biii + -^ran- And for a giveu survey the minimal 
separation ai is determined by the PSF width, whilst the 
maximum separation 02 should be chosen to maximize S/N. 

We estimate the S/N expected from one exposure at 
the Advanced Camera for Surveys (ACS) on board of the 
Hubble Space Telescope (HST). We assume that the sepa- 
rations are distributed according to Opik's law (1924), i.e. 
P{a) oc a~^ , and that the fraction of wide binaries with 
separations 10~^ < (a/lpc) < 1 is 10%, i.e. similar to that 
found in the stellar halo of the Milky Way (Longhitano & 
Binggeli 2010). The stellar luminosity functions were taken 
from Dartmouth Stellar Evolution Program (Dotter et al. 
2008) assuming the Chabrier mass function (Chabrier 2001) 
with Mc = 0.08 Mo and a = 0.952 (Bastian et al. 2010) 
and old (12 Gyr) metal-poor ([Fe/H]=-2) stellar popula- 
tions. Given the PSF width of ACS (Jee et al. 2007), we 
use a minimal separation of 4x0.05" pixels. 

Figure [4] shows the resulting signal-to-noise for dwarf 
galaxies at different heliocentric distances. Diamonds show 
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the expected S/N value for MW dwarf satellites adopting 
miim = 29™. The bottom panel shows the range of binary 
separations probed by HST as a function of distance, where 
ai is set by resolution and 02 is selected to maximize the 
S/N ratio. 

This Figure shows two interesting points. First, the sur- 
face brightness of the dwarf is not the driving parameter for 
the detection of binaries. Instead the most relevant quanti- 
ties are the distance to the dwarf and the depth of the ob- 
servations. Second, the range of binary separations probed 
by HST seems to be reasonable to test the predictions of our 
paper, but in order to put significant limits on the binary 
fraction in dwarfs, several deep exposures are needed. Note 
that Figure HI shows the results of a single ACS exposure. 
For 10 deep exposures the detection signal would boost to 
5.8, 5.0, 4.5, 4.4, 4.2 and 4.0 for Coma, Ursa Minor, Bootes 
I, Ursa Major II, Sculptor and Draco, respectively. 



6 SUMMARY 

The present cosmological paradigm, ACDM, predicts the 
presence of dark matter substructures devoided of baryonic 
material (the so-called "dark galaxies" ) in all galactic haloes 
on all mass scales. We have considered in this paper the pos- 
sibility of using wide stellar binaries to uncover the existence 
of these objects in dwarf spheroidal galaxies (dSphs), whose 
estimated dark matter densities surpass those of any known 
galaxy in the Universe. 

Our analytical estimates suggest that a large fraction of 
wide binaries can be wiped out due to tidal encounters with 
dark substructures. These interactions occur in a "catas- 
trophic" regime, wherein the disruption rate is proportional 
to the density of perturbers averaged over the life-time of 
the binary population (pp). 

We find that in this regime the perturbed binary dis- 
tribution has a separation function that scales as dN/da oc 
a~^'^ at a^Omax, where amax oc {pp)~^'^. Our estimates 
show that for the range of dynamical masses and Galacto- 
centric distances spanned by Milky Way dSphs, all dwarfs 
are expected to have a truncation in the separation function 
at Omax ~ O.lpc within the CDM context. 

We examine in which dSphs observations of large sepa- 
ration binaries are feasible with existing instruments, find- 
ing that the ACS camera on board of the Hubble telescope 
would be able to test the predictions enclosed in this paper 
in dSphs that locate at heliocentric distances ;$ 100 kpc, 
even if the binary fraction amounts only 10% of the stellar 
population. 

Altogether these considerations may pose a strong test 
through which ACDM shall soon have to pass. 
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APPENDIX A: POINT-MASS VS. EXTENDED 
PERTURBERS 

The point-mass approximation may not be appropriate for 
encounters where the binary star penetrates deeply within 
the perturber object, i.e. b/rh^l- On the other hand, in 
!|2] we argue that collisions between dark matter substruc- 
tures and stellar binary systems in the dark matter haloes of 
dSphs are expected to enter within the catastrophic regime. 
In this regime binary stars may be completely shattered by 
single encounters with dark substructures if b^bmin, where 
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Figure Al. Binary separation distribution in simulations where 
dark structures are treated as point masses (upper panel, see 
Fig. ml and extended objects (lower panel). In both panels the 
binary semi-major axis is normalized to the binary maximum 
separation amaxi as computed from eq. ^. Note that the per- 
turbed binary population scales as dNi,/da oc a~^l regardless of 
the inner structure of the perturber objects (dashed line). 



&min is defined in eq. (pi. Thus, fatal encounters with point- 
mass perturbers will occur within an interval of impact pa- 
rameters Vh < b < bmin- 

In the case of large, extended perturbers, however, we 
may find situtations where r^ ^ &min. The goal of this Ap- 
pendix is to examine how the point-mass approximation may 
affect our results. 

Let us consider first how the mass and the size of dark 
structures are related in CDM. Given that the mass and cir- 
cular velocity of CDM haloes scale as M cc V^ , and that 
the virial theorem is, V^ oc M/R, it follows that the half- 
mass radius and the mass of dark substructures scale as 
rn oc mI^^. From eq. ^ we find femin/r,. oc Ml''^a^''^ /rn ~ 
Mp a^'"*, which suggests that whether or not the point- 
mass approximation is valid is mainly set by the binary sep- 
aration. 

Scaling eq. H to the parameters typically found in 
dSphs (see !|2| we find 



I / Mp \/2 ■ lOkms" 



V 108 Me 



Kc 



1/2 



IMq O.lpc 
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1/4 



Thus, for bmin/a ~ 2 ■ 10'' we expect the point-mass 
approximation to hold for binaries with very wide separa- 
tions, a^ (fomin/a) r/i ~ 5- 10~*rh. Conversely, if we adopt 
a fiducial separation a = 0.1 pc, the point-mass approxima- 
tion will be valid for perturbers with half-mass radii Vh S 200 
pc. 

Unfortunately, there is no clear CDM prediction about 
the size of sub^haloes orbiting in dwarf-size haloes, simply 
because no cosmological N-body simulation has yet reached 



Figure A2. Separation at which the binary distribution trun- 
cates as a result of encounters with point-mass (closed symbols) 
and extended (open symbols) substructures that amount a to- 
tal mass of /subAfviri where Mvir is the virial mass of the dwarf 
galaxy. The value of a^g^^ is defined as the separation at which 
dNf,/dlogiQ a drops 20% with respect to the unperturbed distri- 
bution. Dashed lines show the expected position of the truncation 
as computed from eq. ^. The dotted line is a fit to the values of 
"max derived from our N-body simulations. 



the resolution required to study the internal structure of 
these systems in detail. Here we shall assume that the size- 
mass relation found for subhaloes (Diemand et al. 2008; 
Springel et al. 2008) also applies to sub^haloes, and that 
both can be approximated by 



rh 



lO^pc 



M„ 



109 M< 



1/3 



(A2) 



In the upper and lower panels of Fig. |A1| we show the 
results of repeating the simulations presented in !|3] adopting 
point-mass and extended sub^halo models, respectively. In 
practice, the only difference between both calculations cor- 



responds to the value of U(b/rh) adopted in eq. ([T|, which 
U = 1 for point-mass perturbers, and U{x) = a;'^^(l-|-a;)^ 



is 

2.5 



for extended ones (see §8.4 of Binney & Tremaine 2008). 

As already shown in !|3] if dark substructures are as- 
sumed to be point-mass objects eq. ([5| provides an accurate 
(^1) estimate of the binary separation at which deviations from 
the unperturbed binary distribution, dNt/da oc a~^ , start 
to become evident. At a separation a — flmax we find that 
the number of bound binaries drops 20% with respect to 
the unperturbed distribution independently of the perturber 
density. 

Interestingly, if dark substructures are modelled as ex- 
tended objects the truncated binary population also scales 
as dNt/da oc a~^'^ for a 3> Omaxi where we define aJnax ^-s 
the separation at which the number of binaries drops by 
20%. However, in this case eq. (|5| underestimates the lo- 
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cation of the truncation of the binary distribution, so that 

Fig. |A2| shows how the maximum separation, a'^^^, 
varies as a function of the mean sub^halo density (pp) — 
/aubMvir3/(47ri?j) for point-mass (closed symbols) and ex- 
tended (open symbols) bodies. It is clear that the point-mass 
approximation provides an accurate estimate of a'^^y. if the 
fraction of sub^haloes represent a small fraction of the dwarf 
galaxy mass. This is because the number of massive sub- 
structures, which also tend to have larger half-light radii, is 
proportional to /sub- For /sub ~ 1, however, the point-mass 
approximation under-estimates aj^ax by a factor « 10. 

Interestingly, Fig. |A2| also show that a^^^ a-lso scales as 
a power-law if perturbers are modelled as extended objects. 
In particular we find that oj^ax ck {Pp)~^ , a- result that 
is used in !|4] to estimate the truncation in the separation 
function of the binary population orbiting in Milky Way 
dSphs. 
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